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Abstract. Many real life networks present an average path length logarithmic with the number of nodes and 
a degree distribution which follows a power law. Often these networks have also a modular and self-similar 
structure and, in some cases - usually associated with topological restrictions- their clustering is low and 
they are almost planar. In this paper we introduce a family of graphs which share all these properties and 
are defined by two parameters. As their construction is deterministic, we obtain exact analytic expressions 
for relevant properties of the graphs including the degree distribution, degree correlation, diameter, and 
average distance, as a function of the two defining parameters. Thus, the graphs are useful to model some 
complex networks, in particular technological and biological networks. 

PACS. 89. 75. He Networks and genealogical trees - 05.45.Df Fractals - 89.75.Fb Structures and organi- 
zation in complex systems 



1 Introduction 

Ten years have past since the publication of the ground- 
breaking papers by Watts and Strogatz [1] on small-world 
networks and Barabasi and Albert [2] on scale-free net- 
works. Their works led researches to the design of new 
network models to describe complex systems in nature 
and society like the Internet, protein-protein interactions, 
transportation systems or social and economic networks. 
Their models try to match observational studies which 
have identified at least three important common charac- 
teristics for real-life networks: They exhibit a small aver- 
age distance and diameter (compared to a random net- 
work with the same number of nodes and links); the num- 
ber of links attached to the nodes obeys a power-law dis- 
tribution (the networks are scale- free); and recently it has 
been discovered that, often, real networks are self-similar, 
see [3] sometimes showing a degree hierarchy related to 
the modularity of the system. 

Many of the proposed models are stochastic as this 
is the case for the now classical preferential attachment 
method [5]. Thus, the use of probabilistic techniques is re- 
quired to estimate the main parameters of a network [4]. 
However, a deterministic approach has proven useful to 
complement and enhance the probabilistic and simulation 
techniques. Deterministic models have a clear advantage, 
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as they allow an analytical exact determination of rele- 
vant network parameters, which then may be compared 
with experimental data coming from real and simulated 
networks 

Among the different techniques to generate determin- 
istic models are of particular interest those based in recur- 
sive or iterative methods, where at each generation step 
nodes are connected to a known substructure of the net- 
work. This is the case of the pseudo-fractal networks [5] 
where, at each step, all existing links are considered and a 
new vertex is added to each of them. This construction can 
be generalized if cliques of a given size are used instead of 
links (which are 2-cliques) , see [6] . Similar rules give the in- 
teresting Apollonian networks [HElinj . A related technique 
produces networks by duplication of certain substructures, 
see [To). 

A generalization of these former methods introduces at 
each iteration a more complex substructure than a single 
node which is added to the network in a deterministic way. 
Substructures considered are triangles [11] , circles [12] and 
paths [13]. 

In this paper we go an step further in the generaliza- 
tion by considering the introduction of d parallel paths. 
The result is a family of planar, modular, hierarchical and 
self-similar networks, with small-world scale-free charac- 
teristics and with clustering coefhcient zero, and all these 
parameter are determined by d as well as by the iteration 
step t. We note that some important real life networks, 
for example those associated to electronic circuits, Inter- 
net and some biological systems [I4 ] ll5 j . have these char- 
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actcristics as they arc modular, almost planar and with 
a reduced clustering coefficient and have small-world and 
scale-free properties. Thus, these networks can be modeled 
by our construction which can be considered as a new tool 
to study their associated complex systems. In the next sec- 
tion we introduce the family of graphs object of study and 
in Section [3] we calculate analytically some relevant prop- 
erties for the graphs, namely, the degree distribution, de- 
gree correlations, the diameter and the average distance. 
The last section provides some conclusions. 



2 Generation of the graphs Md(t) 

In this section we introduce a family of modular, self- 
similar and planar graphs which have the small- world pro- 
perty and are scale-free. The family depends on an ad- 
justable parameter d and the iteration number t. We pro- 
vide an iterative algorithm, and also a recursive method, 
for its construction. The construction methods allow a di- 
rect determination of the order and size of the graph. 

Iterative construction.- We give here an iterative for- 
mal definition of the proposed family of graphs, Md{t), 
characterized hy t > 0, the number of iterations and a 
parameter d associated with the self-repeating modular 
structure. 

First, we call generating edge the only edge of Md{0) 
and all edges of Md{t) whose end vertices have been in- 
troduced at different iteration steps t. All other edges of 
Md(t) will be known as passive edges. A generating edge 
becomes passive after its use in the construction. 

The graph Md{i) is constructed as follows: 

For t = 0, Md{Q) has two vertices and a generating 
edge connecting them. 

For t > 1, Md{t) is obtained from Md{t—1) by adding, 
to every generating edge in Md{t — 1), d parallel paths P4 
of length three by identifying the two final vertices of each 
path with the endvertices of the generating edge. 

The process is repeated until the desired graph order 
is reached, see Fig. [TJ We note that the graph order can 
be also adjusted with the parameter d (number of parallel 
paths that arc attached to each generating edge). 

Recursive modular construction. - The graph Md(t) is 
also defined as follows: 

For t = 0, Md{0) has two vertices and a generating 
edge connecting them. 

For t = 1, Md{l) is obtained from Md{Q) by adding 
to its only edge d parallel paths P4 of length three by 
identifying the two final vertices of each path with the 
endvertices of the initial edge. 

For t >2, Md{t) is made from 2c? copies of Md{t — 1), 
by identifying, vertex to vertex, the initial edge of each 
Md{t — 1) with the generating edges of Md(l), see Fig. [TJ 

Order and size of Md{t).- We use the following nota- 
tion: V{t), E{t) and Eg{t) denote, respectively, the set of 
vertices, edges and generating edges introduced at step i, 
while V{t) and E(t) denote the set of vertices and edges 
of the graph Md{t). 



'• •' 




Fig. 1. Graphs Md{t) produced at iterations t = Q,1,2 and 3 
for d = 2. 

Notice that, at each iteration, a generating edge is re- 
placed by 2d new generating edges and d passive edges. 
Therefore: \Eg{t + l)\ = 2d - \Eg{t)l and \Eg{t)\ = (2d)*. 
As each generating edge introduces at the next iteration 
2d new vertices and 3d new edges we have 1)| = 2d- 

\Eg{t)\ = (2d)*+i and \E{t + l)\=M-\Eg{t)\ =3d-(2d)*. 
As |y(0)| = 2 and |i?g(0)| = 1, the order and size of M(t), 
t > 0, is: 

i=0 

, , , ~, 3d(2d)* - d - 1 
\E{t)\=Y,\Ei^)\^ . (1) 

Planarity. - A graph is planar if it can be drawn on the 
plane with no edges crossing. By construction of Md(t), 
the introduction at each iteration of d parallel paths con- 
nected to each generating edge, which afterwards becomes 
passive, adds 2d new vertices to the graph and they can 
be drawn without crossing edges. Planarity could also be 
proven from Kuratowski's theorem or from the known pla- 
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narity test which states that a graph is planar if it has no 
cycles of length 3 and \E\ < 2\V\ - 4, \V\ > 3, see [E]. 



3 Topological properties of Md{t) 

Thanks to the deterministic nature of the graphs Md{t), 
we can give exact values for the relevant topological pro- 
perties of this graph family, namely, the degree distribu- 
tion, degree correlations, the diameter and the average 
distance. 

Degree distribution. - Initially, at i = 0, the graph has 
two vertices of degree one. When a new vertex i is added to 
the graph at iteration ti, this vertex has degree 2 and it is 
connected to only one generating edge. We use the follow- 
ing notation: kg{i,t), kp{i,t) and k{i,t) are, respectively, 
the number of generating edges, passive edges and total 
edges connected to vertex i, at step t > ti. Therefore 
k{i,t) = kg{i,t) + kp{i,t) is the degree of vertex i at this 
step. 

From the construction process we can write, 



kp{i, t 4- 1) = kp{i, t) + kg{i, t) 
kg{i, t -I- 1) = dkg{i, t) 



with the initial conditions, 

kg{i,ti) = l ^ti and kp{i,ti) 
we have for > 1. 



if ti^O 

1 otherwise 



kg{i,t) 
kp{i, t) 



■ *' yti and 

PhJ (4) 

2 + otherwise. 



All the vertices that have been introduced at step ti 
have the same degree at step t: 

1. The two vertices introduced at step ti ^ have degree. 



k{i, t) = kg{i, t) + kp{i, t) = 

-d - 1 



+ 1 



1 



1 



(5) 



2. The 11^(^1)1 = (2(i)*' vertices introduced at step ti > 
have degree. 



k{i, t) = kg{i, t) -f kp{i, t) 



jt-ti 



+ 2- 



d-1 



= 1 



d-1 



(6) 



Therefore the degree spectrum of the graph is discrete 
and to relate the exponent of this discrete degree distri- 
bution to the power law exponent of a continuous degree 
distribution for random scale free networks, wc use the 



technique described by Newman in |15j to find the cumu- 
lative degree distribution Pcum(fc)- If we denote by V{t, k) 
the set of vertices that have degree k at step t, 



(2rf)*+i+2d-2 
2d-l 



\v{t)\ 

_ (2d)*'+i + 2d- 2 _ 
~ (2d)*+i + 2d - 2 ^ 

ln(fc + l^-l) 

_ (2d)* +^ + 2d - 2 

~ (2d)*+i + 2d - 2 ■ 

Fot t large, we obtain, 



-Pcum(fc) « (2d) ^ 



(fc 



2-k 



— 1) 



'"(M) 1 2-d ln(2d) 

= k (1 - - + — ;- >°''') . 

d fed 
For fc >> 1 this expression gives 

ln(2d) 1 ln(2d) 

Pcum(fc) «fc-T^(l- -)-T^ 

(2) The degree distribution follows a power-law 

-Pcum(fc) ~ k^'^ with exponent 7 = ^-^^^ see Fig. [2] 



(3) 



(7) 



0.01 




0.001 - 



10-*- 



Fig. 2. Log-log representation of the cumulative degree dis- 
tribution for Mdit), d — 2,3,4,5. The reference line has slope 
-2. 



Therefore the degree distribution is scale- free. Research 
on networks associated to electronic circuits show that 
many of them are almost planar, modular and have a small 
clustering coefficient and in most cases their degree distri- 
butions follow a power-law [141ll5j with exponent values 
in the same range than those of Md{t). 

Correlation coefficient.- We obtain here the Pearson 
correlation coefficient, r{d,t), for the degrees of the end- 
vertices of the edges of Md{t) |17j . 



r(d,0 = 



\m\J:^J^k^-[E^-2iJ^+k^)? 

\m\E^hiJf+kf)-[J:.Kn+k^)? 



(8) 
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where ji, ki are the degrees of the end vertices of the ith 
edge, with i = 1, • • ■ , \E{t)\. 

To find r{d, t) we look at the degree distribution of the 
endvertices of the edges in E{ti) at a given step t;. We 
denote by (j, k) an edge connecting vertices of degrees j 
and k. 

The edges introduced at step ti are: 

1. Edges (2, 2), connecting two vertices introduced at step 
ti > 0. There are (2c?)*'/2 edges (a half of the vertices 
introduced at step ti). Notice that there is one edge 
(1, 1) introduced at ti = 0. 

2. Edges {2,k{i',ti)) connecting vertices of degree two, 
introduced at step ti , with all the vertices i' introduced 
at step til with < i,;' < <i — 1. For each vertex i' there 



are k„ 



edges: 



From the two vertices introduced at ti' = 0, see (O, 
there are 2d*'' edges (2, 

From the (2d)*'' vertices introduced at ti' > 0, see ©, 
there are (2d)*''d*'-*'' edges (2, 1 + '^'^Y-i'"^ )- 

Table [T] displays a summary of these results. 
Using the results of Table [TJ we can find the sums: 



^ jih = (4 d + 16 - 51 + 41 d'' - 8 - 

i 

- 3d'^+d^ + d*+^(40 + 8f-80-2*) + 
+ d*+^(-184 - m + 282 ■ 2') + 
+ d*+'*(306 + 74t - 373 ■ 2') + 
+ d*+^( 



+ d'- 
+ d 



,2i+3 



-236 - 64i + 227-2') + 
.6 + 26t-63-2') + 
-12-4t + 7-2') + 



+ d^''^''(10 + 4t) + d^'+''(-43 - 18t) + 



+ d^'+''(62 + 28t) + d^'+''(-35 - ISt) + 

(2d^ -5d+2)d(d-2)), 
-2 



(-2-3d + 



+ lOd 
+ 

+ d 
- 2d 



(2d- l)(d-2)(d- 1)2 

6d^ + d* + d*+^ (-4 + 16 • 2') + 
*+^(14 - 37 ■ 2*) + d*+^(14 + 26 ■ 2') + 
'+*(4- 

'2t+4N 



5-2' 



^2t+2 _^ 3^2t+3 , 



{ji + ki) = -2(-8d^ - 32d^ + 186d'' - 282d^ + 



lld^° + 



+ 145d^ + 49d'^ - 96d** + 48d'-* 
+ d" + d'+^(-72 + 160 • 2*) + d*+'*(384 - 
- 728 • 2') + d*+''^(-750 + 1252 • 2') + 
+ d*+'^(606 - 882 ■ 2*) + d*+^(-33 - 39 ■ 2*) + 
+ d*+®(-285 + 456 • 2*) + d*+^(201 - 286 ■ 2*) - 
+ d*+^°(-57 + 74 ■ 2*) + d*+"(6 - 7 ■ 2') + 
+ 4d^*+^ - 16d^'+'^ + 17d3'+' + - 



2d=''+")/((d -2d+l) 
(2d - l)(d - 2)(d^ - 2d^ - 2d + 4)d^{d - 1)^). 



Replacing these sums into equation ([S]) we obtain di- 
rectly, and for any d, a (long) exact analytical expression 
for the Pearson correlation coefficient of Md{t). For d = 2 
the equation is displayed in Eq.(9): 

Table [2] shows numerical values of the correlation for 
different instances of the graphs. 





t = 1 


t = 2 


t = 3 


t = 10 


d = 2 


-0.1667 


-0.0886 


-0.0460 


-0.0003 


d= 10 


-0.4091 


-0.2338 


-0.1174 


-0.0009 


d = 100 


-0.4901 


-0.2057 


-0.0934 


-0.0007 



Table 2. Correlation coefficient at steps t — 1,2, 3, 10 for se- 
veral values of d. 



From the values of the correlation coefficient we see 
that this family of graphs has the degrees of the endver- 
tices negatively correlated, large degree vertices tend to 
be connected with low degree vertices, and the graphs are 
disassortative. 

We notice that most technological and biological net- 
works have this property, see [15] . 

For d >> 1, we obtain 



r(d,t) 



-4d4 



1 



d4t+i2(_4 + 6-2*) 
which for t large gives r(d, t) ^ 0. 



1-3-2* 



Diameter.- At each iteration step we introduce, for 
every generating edge, 2d new vertices. These vertices are 
among them at distance at most 3. As each vertex joins 
the graph of the former step through one new edge, the 
diameter will increase by exactly 2 units. Therefore D{t) = 
D(t-l) + 2, t > 2. As ^(l) = 3, we have that the diameter 
of Mrf(t) is D{t) = 3 + 2 • (t - 1), i > 1. Therefore, from 
Eq.[l] and as for t large, t ^ In \V{t)\ we have in this limit 
that D{t) r^ln\V{t)\. 

Average distance.- The average distance of Md(t) is 
defined as: 



1 



\V{t)\i\V{t)\-l)/2 



(10) 



where dij is the distance between vertices i and j. In what 
follows, S{t) will denote the sum J2ijev{t) 
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Step ti 


Edges at step ti 


Number 


Edges at step t > ti 




(1,1) 


1 


\ d-l ' d-1 / 


1<U <t 


(2,2) 

(2,^rr^) 


(2d)'. 
2 

2d'' 


/I 1 dd'^'i-1 1 1 dd*^*i-l\ 
\ d— 1 , d— 1 ' 

/-f 1 dd*"'i-l d'+l-l\ 
V-"- + d-1 ' d-1 / 


2<U <t 
l<ti<U-l 


(2,1 + ^^4^) 


(2d)*- ■ d*'-*'. 


/I 1 dd*^*?— 1 1 1 dd'^'i— 1\ 
N-"- ^ d-1 ' -'- ' d-1 / 



Table 1. Number of edges in Mdit) according to the degrees of their endvertices. 



4*^2 _ 2'+it + 3 ■ 2^*+^ - 2^'+^t + 13 ■ 4* - 3 ■ 2*+^ + 4^'+^ - 3 ■ 2^'+^ + 1 
24t+if2 22'+it - 23'+:'f + 2**+3f - 2« + 5 ■ 4« - 23*+'i + 3 ■ 24'+3 - 3 . 25«+2 




copies of Mdit). 

The modular recursive construction of Md{t) allows us 
to calculate the exact value of D{t). At step t, Md{t + 1) 
is obtained from the juxtaposition of 2d copies of Md{t), 

which we label M^^^ , rj = 1,2, ■■■ ,2d, see Figures [T] and [31 
Whenever possible, we drop the subscript d and represent 
I as Ml to keep the notation uncluttered. The copies 
are connected one to another at 2d + 2 vertices which 
we call connecting vertices such as u, v, w, x, y, and z 
in Fig. [31 Thus, the sum of distances St+i satisfies the 
following recursion: 

St+i=2dSt + At. (11) 

where At is the sum over all shortest path length whose 
endpoints are not in the same Af/''^ branch. 



To compute At, we classify the vertices of Md{t + 1) 
into two categories: the two vertices with the largest de- 
gree (i.e., w and x in Fig. [3]) are called hubs, while any 
othe vertex is named a non-hub vertex. Thus At can be 
obtained by summing the following path lengths that are 
not included in the distance between vertex pairs ofM.^"): 
length of the shortest paths between non-hub vertices, 
length of the shortest paths between a hub and non-hub 
vertices, and length of the shortest paths between hubs 
(for example, dux, and duz)- 

Let us denote A^''^ as the sum of all shortest paths be- 
tween non-hub vertices, whose endpoints are in A^"^ and 
Mt^\ respectively. Thus, A'^'^ rules out the paths with 
endpoints at the connecting vertices belonging to or 
Mt^\ For example, each path contributing to A\''^ does 
not end at vertex u, v, w or x, and each path contributing 
to A]''^ does not end at vertex u, w, x or z. According to 
its value, A"''^ can be split into three classes, where the 
three representatives are AI'^ , A]'^, and A]''^, and the car- 
dinality of the three classes are d, d{d — 1), and d{d — 1), 
respectively. Analogously, the length of the shortest paths 
between a hub and all non-hub vertices can be classified 
into two classes, while the shortest paths between hubs 
can be partitioned into three classes with path lengths 
equal to 1, 2, or 3. 

Let be the set of non-hub vertices in Mt"\ then 
the total sum At is given by 

At = dA]^^ + d{d - 1) (^Al^^ + A]''^^ + 2d{d + 1) 

d.u,j + 2d{d - 1) 51 + 

+ {d + l)duv + d{d + 1) duy + d{d - 1) duz, (12) 
where duv ~ 1, duy ~ 2, and duz ~ 3 arc easily seen. 
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Having At in terms of the quantities of A]'^. A]'^, 
gj^2 dwj J2j£n^ duj, and J^j^n^ the next step 
is to explicitly determine these quantities. To this end, we 
classify non-hub vertices in Md{t +1) into two different 
parts according to their shortest path lengths to either 
of the two hubs (i.e. w and a;). Notice that the vertices w 
and X themselves are not partitioned into either of the two 
parts represented as Pi and P2, respectively. The classi- 
fication of vertices is shown in Fig. [3]). For any non-hub 
vertex (p, we denote the shortest path length from ip to 
w, X as a, and 6, respectively. By construction, a and h 
can differ at most by 1 since vertices w and x are adja- 
cent. Then the classification function class(ip) of vertex ip 
is defined to be 



class(ip) 



Pi 
P2 



for a < b, 
for a > b. 



(13) 



It should be mentioned that the definition of the vertex 
classification is recursive. For instance, class Pi and P2 in 

Mj^-* belong to class Pi in Md{t + 1), class Pi and P2 in 
(2) 

Ml belong to class P2 in Md{t-\- 1), and so on. Since the 
two hubs w and x are symmetrical, in the graph we have 
the following equivalent relations from the viewpoint of 
class cardinality: classes Pi and P2 are equivalent one to 
each other. We denote the number of vertices in network 
Md{t) that belong to class Pi as Nt^p^, and the number 
of vertices in class P2 as Nt^p^- By symmetry, we have 
Nt,Pi = -^t,P2: which will be abbreviated as Nt hereafter. 
It is easy to see that 



Nt 



\v{t)\ 



1 



d{2dY 



2d - 1 



(14) 



For a vertex ip in Md{t -1-1), we are also interested in 
the smallest value of the shortest path length from p to 
either of the two hubs w and x. We denote the shortest 
distance as this value by and it can be defined as 



= min{a, b). 



(15) 



Let 5t.p^ {^t,P2) denote the sum of for all vertices 
belonging to class Pi (P2) in Md{t). Again by symmetry, 
we have St^p^ = ^t.Pa that will be written as St for short. 
Taking into account the recursive method of constructing 
Md{t), we notice that the vertex classification follows also 
a recursion. Therefore we can write the following recursive 
formula for St+i'. 



St+i=2d5t + dNt + d. 



(16) 



Substituting equation (|14p into equation (jl6[) , and consid- 
ering the initial condition So = 0, equation (flG]) is solved 
inductively 

_2d-2d^ - (2d)i+* + d(2d)i+* - dt{2dY + dt{2dY+* 
' ~ 2(2d- 1)2 ■ 

(17) 

We now return to compute equation (|12p . For conve- 
nience, we use Ft'^ to denote the set of non-hub vertices 



belonging to class Pi in m/''\ Then Z\('^ can be written 



as 



Z\^2 = ^ 4 



^ ^ (j^rw dyjx 4" dxs) 4~ ^ ^ {drw dyjy -h dvs) 
{dru + dux +dxs) + ^ {dru + duv + d.us) 

(18) 



serf 



serf 



= 8NtSt + 6iNt)^. 
Analogously, we find 

and 



(19) 



(20) 



Al''' = 8NtSt + 8iNty. 

Next we will determine other quantities in equation ([12 
with T^jen^ '^wj given by 



^2St + 3Nt. 
Analogously, we can obtain 



(21) 



(22) 



duj =2St + bNt. 



Substituting equations jSO]), ^ and ^ 

into equation (|12p . we have the final expression for cross 
distances At, 

At^l + U^ + 4rf(4d - l)Nt + Qd{2d - l){Ntf + 
+ MSt[d+{2d-l)Nt] = 

^ -2d^ + {l-d){2df+^ + 



(1 - 2kf 

+ + 2i)4i+*d4+2* - (7 + 2t)d^{2df+'^'). 



(23) 



Inserting equation (j23p into equation (jlip and using 
the initial condition = 1, equation pT|) is solved induc- 
tively. 



St = 



1 



r(-l + Ad- + 2d^ + 2i+*di+* 



(-l + 2d)3' 
- 7 • 2^*^2+2* + 3 • 2^+2*^3+2* _ 2i+*di+*t + 
+ 3 • 2^+'d^+h - 2^+'d'+H - 21+2*^2+2*^ + 

+ 22+2*d3+2tt). (24) 

Substituting equation into equation (|10p yields the 
exact analytic expression for the average path length of 
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Md{t) as 

D{t) = (-1 + Ad- + 2d^ + 2^+*d^+* - 7 ■ 2^*^2+2* _^ 
+ 3 • 2i+2*d3+2* - 2'+'d'+H + 3 • 2'+'d^+H - 

_ 22+trf3+tt _ 2^+2*^2+2*^ + 22+2*^3+2*^) 

/ ((-1 + 2d){-l + d + 2*di+*)(-l + 2i+*di+*)). 

(25) 

Notice that for a large iteration step, t — > oo, D{t) ~ 
t ~ ln|T^(f)|, which shows a logarithmic scaling of the 
average distance with the order of the graph. As wc have 
a similar behavior for the diameter, the graph is small- 
world. 



4 Conclusion 

The graphs Md{t) introduced and studied here are planar, 
modular, have a disassortative degree hierarchy and are 
small-world and scale-free. Another relevant characteris- 
tic of the graphs is their clustering zero. A combination of 
a low clustering coefficient, modularity, and small-world 
scale-free properties can be found in some real networks, 
in particular in technical and biological networks [T51IT4]. 
and most of them are also disassortative. As examples, 
the largest benchmark considered in |14| -a network with 
24097 nodes, 53248 edges, average degree 4.34 and average 
distance 11.05- has a degree distribution which follows a 
power-law with exponent 3.0, and it has a small clustering 
coefficient C = 0.01 and other network properties (diame- 
ter, average distance, average degree, etc.) are also in the 
same range than those of the graph Mq{A), see [15] . Also, 
the S. cerevisiae protein-protein interaction network has 
1870 and 2240 edges, see [18], and again most of its net- 
work parameters, as described in [15] can be compared 
directly with those of the graph Me{3). 

Finally, wc should emphasize that the planar property 
and the deterministic character of the family, in contrast 
with more usual probabilistic approaches, should facili- 
tate the exact determination of other network parameters 
and the development of new network algorithms that then 
might be extended to real-life complex systems. 
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